Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 24 November 2010 



(MN I^TbX style file v2.2) 



Non-stationary dynamo and magnetospheric accretion 
processes of the classical T Tauri star V2129 Oph 

o ■ 

g : J.-F. Donati 1 *, J. Bouvier 2 , F.M. Walter 3 , S.G. Gregory 4 , M.B. Skelly 1 , 

C\i : G.A.J. Hussain 5 , E. Flaccomio 6 , C. Argiroffi 7 , K.N. Grankin 8 , M.M. Jardine 9 , 

> 'F. Menard 2 , C. Dougados 2 , M.M. Romanova 10 & the MaPP collaboration 

O ' 



m 

(N 

Oh 
i 

O 



(N 
> 
OS 

oo 



1 LATT-UMR 5572, CNRS & Univ. de Toulouse, 14 Av. E. Belin, F-31400 Toulouse, France 

2 LAOG-UMR 5571, CNRS & Univ. J. Fourier, 414 rue de la Piscine, F-S8O4I Grenoble, France 

3 Department of Physics and Astronomy, Stony Brook University, Stony Brook NY 1 1794-3800, USA 

4 School of Physics, Univ. of Exeter, Stacker Road, Exeter EX4 4QL, UK 

5 ESO, Karl-Schwarzschild-Str. 2, D-85748 Garching, Germany 

6 INAF, Osservatorio Astronomico di Palermo, Piazza del Parlamento 1, 90134 Palermo, Italy 

7 DSFA, Univ. di Palermo, Piazza del Parlamento 1, 90134 Palermo, Italy 

8 Crimean Astrophysical Observatory, Nauchny, Crimea 334413, Ukraine 

9 School of Physics and Astronomy, Univ. of St Andrews, St Andrews, Scotland KY16 9SS, UK 

10 Department of Astronomy, Cornell University, Ithaca, NY 14853-6801, USA 



2010 November, MNRAS in press 



ABSTRACT 

We report here the first results of a multi- wavelength campaign focussing on mag- 
netospheric accretion processes of the classical T Tauri star (cTTS) V2129 Oph. In 
this paper, we present spectropolarimctric observations collected in 2009 July with 
ESPaDOnS at the Canada-France-Hawaii Telescope (CFHT) and contemporaneous 
photometry secured with the SMARTS facility. Circularly polarised Zeeman signatures 
are clearly detected, both in photospheric absorption and accretion-powered emission 
lines, from time-series of which we reconstruct new maps of the magnetic field, photo- 
spheric brightness and accretion-powered emission at the surface of V2129 Oph using 
our newest tomographic imaging tool - to be compared with those derived from our 
old 2005 June data set, reanalyzed in the exact same way. 

We find that in 2009 July, V2129 Oph hosts octupolar and dipolar field compo- 
nents of about 2.1 and 0.9 kG respectively, both tilted by about 20° with respect to 
the rotation axis; we conclude that the large-scale magnetic topology changed signifi- 
cantly since 2005 June (when the octupole and dipole components were about 1.5 and 
3 times weaker respectively), demonstrating that the field of V2129 Oph is generated 
by a non-stationary dynamo. We also show that V2129 Oph features a dark photo- 
spheric spot and a localised area of accretion-powered emission, both close to the main 
surface magnetic region (hosting fields of up to about 4 kG in 2009 July). We finally 
obtain that the surface shear of V2129 Oph is about half as strong as solar. 

From the fluxes of accretion-powered emission lines, we estimate that the observed 
average logarithmic accretion rate (in M© yr _1 ) at the surface of V2129 Oph is — 9.2± 
0.3 at both epochs, peaking at —9.0 at magnetic maximum. It implies in particular 
that the radius at which the magnetic field of V2129 Oph truncates the inner accretion 
disc is 0.93 x and 0.50 x the corotation radius (where the Keplerian period equals the 
stellar rotation period) in 2009 July and 2005 June respectively. 

Key words: stars: magnetic fields - stars: formation - stars: imaging - stars: rotation 
- stars: individual: V2129 Oph - techniques: spectropolarimetry 



1 INTRODUCTION 



* E-mail: donati@ast.obs-mip.fr 



Magnetic fields are thought to play an important role dur- 
ing the formation of stars and their planetary systems. In 
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particular, they presumably oppose the collapse of the gi- 
ant molecular clouds from which stars and planets form, 
strongly inhibit its fragmentation into multiple globules, ef- 
ficiently evacuate most of the accreting mass (e.g., through 
conical outflows and collimated jets) and rapi dly dissipate 
the in itial reservoir of angular momentum fe.g-. lAndre et abl 
2008); they also likely affect the formation and migration of 
planets in the accretion discs surrounding newly-born stars. 

The impact of magnetic fields on the formation of low- 
mass stars (i.e., with M* < 3 Mq) that are still surrounded 
with a gaseous accretion disc - the so-called classical T Tauri 
stars or cTTSs - has been the subject of intense scrutiny and 
vivid debates in the last two decades. With their intense 
large-scale fields, cTTSs are putatively capable of disrupt- 
ing the inner regions of their accretion discs and of braking 
their rotation rather drastically, although the exact physi- 
cal process through which they can re ach this result is still 
controversial (e.g.. iBouvier et aL I l2007l . for a review). 

Following several papers reporting the p resence of 
strong fields on cTTSs (e.g., Ijohns-Krull I2007T I, the first 
spectropolarimetric studies based on time-series of Zee- 
man signatures from accreting and non-accreting regions 
at the surfaces of cTTSs revealed that low-mass proto- 
stars indeed possess large-scale fields, but that the inten- 
sity and topologies of their fields strongly depend on the 
protostellar mass in a way ve ry reminiscent to those of 
low-mass main-sequence stars dDonati et al ] 120071 . SpOSa; 
iDonati fc Landstreed[2009l ; iDonati et al.ll2010al '). These ini- 
tial results however call for a confirmation through surveys 
on a statistically significant sample of cTTSs, allowing in 
particular to study how the properties of their large-scale 
magnetic fields depend on stellar parameters such as mass, 
age, rotation and accretion rates. Ultimately, such investiga- 
tions can guide us towards more realistic models of forming 
Suns, taking into account the effect of magnetic fields in 
particular. 

Magnetic Protostars and Planets (MaPP) is an inter- 
national project focussing specifically on this issue. It has 
been granted 690 hr of observing time over 9 consecu- 
tive semesters (2008b to 2012b) with the ESPaDOnS spec- 
tropolarimeter on the 3.8-m Canada-France-Hawaii Tele- 
scope (CFHT) to survey 15 cTTSs and 3 protostellar ac- 
cretion discs of FU Ori type (FUOrs); it also regularly 
benefits from contemporaneous observations with the NAR- 
VAL spectropolarimeter on the 2-m Telescope Bernard Lyot 
(TBL) as well as photometric observations from Crimea, 
Uzbekistan and Armenia. Additional multiwavelength ob- 
servations from space and/or from the ground are also or- 
ganised in conjunction with MaPP campaigns on a few spe- 
cific targets, providing deeper insights into the physical pro- 
cesses under scrutiny (and in particular magnetospheric ac- 
cretion) . 

Following a first MaPP paper concentrating on the pro- 
totypical low-mass cTTS AA Tau and suggesting that its 
large-scale magnetic field is strong enough to expel out- 
wards (e.g., through a conical wind) most of the accreted 
disc material (the prope ller regime) and thu s to spin-down 
the protostar efficiently (|Donati et al.|[2010lj ). we present a 
new study revisiting the properties of the large-scale mag- 
netic field of V212 9 Oph (after the intial investigation of 
IDonati et all 120071 ); in addition to being based on a much 
better and longer spectropolarimetric data set, this study 



also benefits from additional contemporaneous multiwave- 
length observations from both space (at X-ray frequen- 
cies with Chandra/HETGS) and from the ground (high- 
resolution spectroscopy from HARPS@ESO, low resolution 
spectroscopy and photometry from the SMARTS facility). A 
brief overview of th e goals of the multiwa velength observing 
program is given in iGregorv et al.l l|2009l ) . 

The present paper concentrates mostly on the spec- 
tropolarimetric data, while the other related data sets (and 
in particular the Chandra/HETGS spectra, Argiroffi et al, 
in preparation) will be described and analysed in forthcom- 
ing companion papers. After briefly recalling the stellar pa- 
rameters of V2129 Oph, (see Sec. [5}, we describe the ob- 
served variations of photospheric lines and accretion prox- 
ies (Sec. [3J and detail their subsequent modelling with our 
magnetic imaging code (Sec.[S]); we also include a complete 
remodelling of our older data to ease the comparison of the 
magnetic images at both epochs. We finally discuss the im- 
plications of these new results for our understanding of mag- 
netospheric accretion processes in cTTSs (Sec. [6}. 



2 V2129 OPH = SR 9 = ROX 29 

V2129 Oph is the brightest cTTS in the p Oph star for- 
mation cloud, with limited visual extinction (~0.6) and a 
photospheric temp erature of ~4500 K (spectral type K5, 
IDonati et alj|2007h. It exhibits strong Ha em ission typical 



to cTTSs (e.g.. IBouvier fc Appenzellerlll992l ) and clear in- 
frared excesses and 10 /im sili cate emission line fl uxes typical 
of dusty accretion discs (e.g., Furlan et al.ll2009l ) suggesting 
a disc mass of ~1 Mq^ l|Cieza et al.ll201ofT 

The distance to p Oph has recently been redeter- 
mined with impr oved accuracy using VLB A, and estimated 
to 120 ± 5 pc (lLoinard et all 120081 '). From the observed 
maximum visual b rightness of V2129 Oph (mv ~ 11.2, 
iGrankin et al.ll2008f). the estimated unspotted visual magni- 
tude (Tl.0±0.2, lDonati et"aill2007h and the vi sual bolomet- 
ric co rrection adequate to V2129 Oph (_ — 1.2. lBessell et alj 
[1998J), one can derive a logarithmic luminosity (with re- 
spect to the Sun) of 0.15 ± 0.1 and thus a radius of R* — 
2.0 ± 0.3 Rq. This is slightly smaller than that derived in 
the previous study, directly reflecting the shorter distance 
to p Oph used here. 

The rotation period of V2129 Oph estimated from reg- 
ular photometric variations is found to be 6.53 d in aver- 
age, varying from about 6.35 to 6.60 d depending on the 
observing season (|Grankin et al.l l2008); when compared to 
other similar stars of the same sample (e.g., V410 Tau, ex- 
hibiting no such behaviour), we conclude that this period 
fluctuation is real. It suggests that differential rotation is 
shearing the spotted photosphere, producing the observed 
period fluctuations as starspots migrate in latitude on long 
timescales - a suggestion that we independently confirm in 
our study (see Sec.[5)l . From the rotational broadening of line 
profiles, the line-of-sight projected equatorial rotation veloc- 
ity of V2129 Oph i s found to be vsini = 14.5 ± 0.3 kms" 1 
l|Donati et al .1120071 ) . where i is the inclination of the rotation 
axis to the line of sight. It implies that i?*sini ~ 1.8 R©, 
assuming that the rotation period at the equator is close 
to the lowest observed periods (i.e., that the photospheric 
shear of V2129 Oph is Sun-like, with an equator rotating 
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faster than the poles, as in most other low-mass stars show- 
ing differential rotation). 

Coupling with the above radius estimate, it suggests 
that i ~ 65°, with potential values ranging from 50° to 
90°. Tomographic imaging (including t his study) suggests a 
smaller inclination angle of i ~ 45 ±15° (|Donati et al.ll2007l) . 
We therefore use an intermediate value, i ~ 60°, slightly 
larger than in the previous study as a result of the shorter 
distance assumed for p Oph. The corresponding value of the 
stellar radius is i?* ~ 2.1 R© which, when coupled to the 
temperature, yields a mass of M* ~ 1 .35 Mq and an age 
of about 2 — 3 Myr (jSiess et aUl200C| ): models suggest in 
addition that V2129 Oph is no longer fully convective and 
hosts a core whose fractional mass and radius are about 
one third those of the star. In this context, the corotation 
radius (at which the Keplerian period is equal to the rotation 
period of the star) is equal to r cor ~ 7.7 R+ or 0.076 AU. 

Mass accretion on V2129 Oph is only moderate. From 
the equivalent widths and the corresponding line fluxes of 
the emission lines usually considered as accretion proxies 
(and in particular the He I D3 line and the Ca n infrared 
triplet, IRT) and using the publis hed empirical cor relations 
between lines and accretion fluxes (Fang et al. 2009]), we can 
derive an estimate of the average logarithmic mass accre- 
tion rate (in M0yr -1 ) at the surface of V2129 Oph, that 
we find to be equal to —9.2 ± O.jQ. Our data also suggest 
that the mass accretion rate of V2129 Oph is more stable on 
the long-term (see Sec. QJ than that of cTTSs like AA Tau 
(showing dominant intrinsic variability on short timescales, 
iDonati et al.l [20 10b); this preliminary conclusion however 
needs to be confirmed with longer data sets, to check in 
particular that we did not catch V2129 Oph in an unusually 
quiet state of accretion. Optical veiling, i.e., the apparent 
weakening of the photospheric spectrum (presumably caused 
by accretion) is apparently weak for V2129 Oph and does 
not exceed a few %. 



3 OBSERVATIONS 

Spectropolarimetric observations of V2129 Oph were col- 
lected from 2009 July 01 to July 14 using ESPaDOnS on 
the CFHT. ESPaDOnS collects stellar spectra spanning the 
whole optical domain (from 370 to 1,000 nm) at a resolving 
power of 65,000 (i.e., 4.6 kms -1 ) and with a spectral sam- 
pling^ 2.6 k m s - , in either circular or linear polarisation 
l|Donati|[200l ). A total of 23 circular polarisation spectra 
were collected, at a rate of 2 spectra per night during the 
first 10 nights. All polarisation spectra consist of 4 individ- 
ual subexposures lasting each 690 s and taken in different 
polarimeter configurations to allow the removal of all spuri- 
ous polarisation signatures at first order. 

All raw frames are processed with Libre ESpRIT, 



1 Note that this new estimate is an order o f magnitude smalle r 
than that mentioned in our previous study llDonati et al.|[200^ . 
at the time simply derived by averaging measurements obtained 
with different methods and thus likely far less accurate; in partic- 
ular, our new mass accretion rate, computed in ex actly the same 
way as those of V22 47 Oph iDonati et al.ll2010aT) and AA Tau 
llDonati et al.ll2010bri , is at least much more consistent with our 
latter MaPP studies. 



Table 1 . Journal of observations collected in 2009 July. Columns 
1 — 4 respectively list the UT date, the heliocentric Julian date 
and UT time (both at mid-exposure), and the peak signal to 
noise ratio (per 2.6 kms" 1 velocity bin) of each observation (i.e., 
each sequence of 4 X 690 s subexposures) . Column 5 lists the rms 
noise level (relative to the unpolarized continuum level I c and 
per 1.8 kms -1 velocity bin) in the circular polarization profile 
produced by Least-Squares Deconvolution (LSD), while column 6 
indicates the rotational cycle associated with each exposure (using 
the ephemeris given by Eq. [TJ . 



Date 


HJD 


UT 


S/N 


CLSD 


Cycle 




(2,455,000+) 


(h:m:s) 




(10 -4 /c) 


(225+) 


T11I 01 

J Hi UT 


-LO. / UOUZ 




1 RO 
loU 


9 r » 
. 


n 

u.uyz 


Jul 01 


13.90207 


09:33:18 


180 


2.8 


0.712 


Jul 02 


14.77090 


06:24:30 


200 


2.2 


0.845 


Jul 02 


14.90375 


09:35:50 


200 


2.2 


0.866 


Jul 03 


15.76152 


06:11:05 


190 


2.3 


0.997 


Jul 03 


15.89192 


09:18:52 


210 


1.9 


1.017 


Jul 04 


16.75643 


06:03:51 


200 


2.1 


1.150 


Jul 04 


16.88363 


09:07:02 


210 


1.9 


1.169 


Jul 05 


17.75614 


06:03:32 


190 


2.3 


1.303 


Jul 05 


17.88444 


09:08:18 


210 


2.0 


1.322 


Jul 06 


18.75731 


06:05:19 


100 


4.4 


1.456 


Jul 06 


18.88404 


09:07:49 


190 


2.2 


1.475 


Jul 07 


19.75701 


06:04:59 


120 


4.2 


1.609 


Jul 07 


19.88577 


09:10:24 


160 


3.1 


1.629 


Jul 08 


20.75692 


06:04:57 


170 


2.6 


1.762 


Jul 08 


20.88477 


09:09:04 


140 


3.4 


1.782 


Jul 09 


21.75541 


06:02:52 


150 


3.0 


1.915 


Jul 09 


21.88237 


09:05:42 


160 


2.7 


1.935 


Jul 10 


22.75503 


06:02:26 


170 


2.4 


2.068 


Jul 10 


22.88263 


09:06:11 


180 


2.3 


2.088 


Jul 12 


24.78050 


06:39:19 


150 


3.0 


2.378 


Jul 13 


25.84672 


08:14:48 


180 


2.7 


2.542 


Jul 14 


26.84682 


08:15:03 


150 


3.0 


2.695 



a fully automatic reduction package/pipeline available at 
CFHT. It automatically performs optimal extraction of ES- 
PaDOnS unpolarized (Stokes I) and circularly polarized 
(S tokes V) spectra gro ssly following the procedure described 
in IDonati et al.l l|l997l ). The velocity step corresponding to 
CCD pixels is about 2.6 kms -1 ; however, thanks to the fact 
that the spectrograph slit is tilted with respect to the CCD 
lines, spectra corresponding to different CCD columns across 
each order feature a different pixel sampling. Libre ES- 
pRIT uses this opportunity to carry out optimal extraction 
of each spectrum on a sampling grid denser than the origi- 
nal CCD sampling, with a spectral velocity step set to about 
0.7 CCD pixel (i.e. 1.8 kms -1 ). All spectra are automati- 
cally corrected of spectral shifts resulting from instrumental 
effects (e.g., mechanical flexures, temperature or pressure 
variations) using telluric lines as a reference. Though not 
perfect, this procedure provides spectra with a relative ra- 
dial velocity (RV) pr ecision of better than 0.030 kms -1 (e.g., 
IDonati et al]|2008bl ). 

The peak signal-to-noise ratios (S/N, per 2.6 kms -1 
velocity bin) achieved on the collected spectra (i.e., the se- 
quence of 4 subexposures) range between 100 and 210 de- 
pending on weather/seeing conditions, with a median of 180. 
Rotational cycles E are computed from heliocentric J ulian 
dates according to the ephemeris of lDonati et alj (|2007t ) : 



HJD = 2453540.0 + 6.53E 



(1) 
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The full journal of observations is presented in Table [T] Ob- 
servations collected on July 04 (rotation cycles 1.150 and 
1.169) were recorded with the (full) moon at an angular dis- 
tance of only about 2° from V2129 Oph. As a result, the 
corresponding unpolarized spectra show significant contam- 
ination and were discarded from the study; circularly po- 
larised spectra are however unaffected and were kept in the 
analysis (see Sec.[4]and Fig. [1]). 

Least-Squares Deconvolution (XSD. lDonati et al.lll997T ) 
was applied to all observations. The line list we employed 
for LS D is computed from an Atlas9 LTE model atmo- 
sphere (|Kuruczlll993h and corresponds to a K5IV spectral 
type (Teff = 4, 500 K and log g = 3.5) appropriate for 
V2129 Oph. Only moderate to strong atomic spectral lines 
(with line-to-continuum core depressions larger than 40% 
prior to all non-thermal broadening) are included in this list; 
spectral regions with strong lines mostly formed outside the 
photosphere (e.g., Balmer, He, Can H, K and IRT lines) 
and/or heavily crowded with telluric lines were discarded. 
Altogether, about 8,500 spectral features (with about 40% 
from Fe i) are used in this process. Expressed in units of the 
unpolarized continuum level 7 C , the average noise levels of 
the resulting Stokes V LSD signatures are ranging from 1.9 
to 4.4xl0 -4 per 1.8 kms -1 velocity bin. 

Almost simultaneous photometry was collected with the 
SMARTS facility in BVRIJHK bands from rotation cycles 
224.3 to 226.9, with a total of 10 to 13 data points depend- 
ing on the band. These data will be analysed in detail in a 
forthcoming companion paper; we only use here the V band 
measurements as a guideline for checking the consistency of 
our modelling. 



4 VARIATIONS OF PHOTOSPHERIC LINES 
AND ACCRETION PROXIES 

As in our previous run l|Donati et al.ll2007t ). Zeeman signa- 
tures are detected at all times in Stokes V LSD profiles, with 
peak-to-peak amplitudes varying from 0.2 to 1.0% (see Fig.[T] 
top right panel). Temporal variations are dominant, the av- 
erage Stokes V signature being most of the time smaller than 
the individual ones; moreover, Zeeman signatures are com- 
plex and feature several reversals throughout the line profile, 
suggesting that the parent field topology is not simple. As 
a result, the line-of-sight projected component of the field 
averaged over the visible stellar hemisphere and weighted by 
brightness inhomogeneities (called longitudinal field and es- 
timated from the first moment of the Stokes V profile, e.g., 
iDonati et al.|[l997l ). varying from —360 to +90 G, is showing 
clear non-sinusoidal fluctuations (see Fig. [2] lower left panel) 
similar to (and about 50% larger than) those observed in 
2005 June. 

LSD Stokes / profiles of photospheric lines, centred at 
an average RV of about — 7 kms , exhibit temporal shape 
variations, and at times clear asymmetries, (e.g., at rotation 
cycles 1.303 and 2.542) very reminiscent to those caused 
by cool spots at the surface of active stars. This is again 
consistent with the conclusions reached from our first set 
of observations, where line profiles distortions (smaller than 
but similar to those reported here) and their evolution with 
time were successfully int erpreted in terms o f cool spots at 
the surface of V2129 Oph |Donati et alj|2007f l. The observed 



photometric variations (see Fig. [2] upper panel) and in par- 
ticular the corresponding colors (not shown here, to be pub- 
lished in a forthcoming companion paper) also strongly ar- 
gue that cool spots are present at the surface of the star, 
with maximum spottedness occurring around phase 0.55, 
i.e., very close to the phase of strongest photospheric longi- 
tudinal field. 

We also find that Stokes I LSD profiles show equivalent 
widths nearly equal to those in the previous data set and 
roughly constant with time (within ~ 2%), suggesting that 
veiling is small. Assuming that the LSD Stokes / profiles 
with maximum equivalent width correspond to an unveiled 
spectrum, we derive that veiling is always smaller than 5% 
(see Fig. [2]), except on July 05 where it reaches about 7% 
as a likely result of a spectral pollution from the nearby 
full Moon. Comparison with spectral lines of non-accreting 
TTSs of similar spectral types (e.g., V410 Tau) confirms this 
conclusion. 

We find in particular that LSD Stokes I and V profiles 
repeat well from one rotation cycle to the next, suggesting 
that rotational modulation clearly dominates over intrinsic 
variability. This is particularly clear on the longitudinal field 
curve (Fig. [2] lower left panel) that repeats almost identi- 
cally over the 2 successive rotation cycles. In particular, this 
behaviour is radically different from that of AA Tau where 
intrinsic variability st rongly dominates the photospheric lon- 
gitudinal field curve (|Donati et al.ll2010bl V 

Emission in the He I D3 line at 587.562 nm, presumably 
produced in the post-shock region at the footpoints of ac- 
cretion funnels, is considered as the most reliable accretion 
proxy. On V2129 Oph, He I emission shows up as a narrow 
profile and amounts to an average equivalent width of about 
11 kms -1 (i.e., 0.022 nm), with values ranging from 8 to 
19 kms -1 (see Fig.[2]and top left panel of Fig. [2]). Maximum 
He 1 emission is found to occur at phase 0.55, i.e., roughly co- 
incident with times of minimum brightness (i.e., maximum 
spottedness); it suggests in particular that the accretion re- 
gion traced by He I emission is close to the cool spot at 
photospheric level. Centred at an average RV of kms -1 , 
He 1 emission is significantly red-shifted (by ~ 7 kms -1 ) 
with respect to LSD photospheric profiles, indicating that 
it forms in a non-static region above the surface (e.g., the 
chromosphere) where the post-shock plasma is decelerating 
towards the star; the line also exhibits sinusoidal RV varia- 
tions of about ±4 kms -1 about the mean, suggesting that 
the corresponding emission concentrates in a region centred 
at phase 0.55 (consistent with the observed maximum He I 
emission) and presumably at high latitudes (to account for 
the small amplitude of the observed RV variations) . At max- 
imum emission, the emission flux is found to vary by almost 
a factor of 2 (from 11 to 19 kms -1 , i.e., 0.022 to 0.037 nm) 
between successive rotation cycles (see Fig. [2}; the disper- 
sion is smaller at emission minimum (around 9 kms -1 , i.e., 
0.018 nm). During our previous run, He I emission fluxes 
were grossly similar, varying from 6 to 19 kms -1 (i.e., 0.012 
to 0.037 nm) across the rotation cycle. 

Clear Zeeman signatures are detected in conjunction 
with He 1 D3 emission (see Fig. [3] top right panel), tracing 
longitudinal fields ranging from 1.5 to 4.5 kG with mag- 
netic maximum occurring at phase 0.55, i.e., at times of 
minimum brightness and maximum He I emission (see Fig. [2] 
second panel left column). As usual for the He I line (e.g., 
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Figure 1. Variation of the Stokes I (left) and Stokes V (right) LSD profiles of the photospheric lines (top) and Ca II IRT lines (bottom) 
of V2129 Oph in 2009 July. The Stokes / profiles of Ca II emission is shown before subtracting the underlying (much wider) photospheric 
absorption, hence the reduced flux in the far wings (with respect to the unit continuum). To emphasize variability, the average profile 
over the run is shown in red. Rotation cycles (as listed in Table 1) and 3cr error bars (for Stokes V data only) are also included next to 
each profile. Unpolarized profiles recorded on rotation cycles 1.150 and 1.169 are significantly affected by the moon and were discarded 
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Figure 2. Temporal variations of the integrated brightness (top panels), He I D3 (second panels), Ca II IRT (third panels) and photo- 
spheric lines (bottom panels) of V2129 Oph in 2009 July. Longitudinal field variations are shown on the left side (3 lower panels) while 
equivalent width and veiling variations are shown on the right side. Data collected within cycles 0.6 — 1.6, 1.6 — 2.6 and 2.6 — 3.6 are 
respectively plotted as red circles, green squares and blue triangles; photometric points shown as black circles were collected one rotation 
cycle before. Fits with sine/cosine waves (plus first overtones for the left panels) are included (and shown as dashed lines) to outline the 
amount of variability caused by rotational modulation. ±1 cr error bars on data points are also shown whenever larger than symbols. 



iDonati et alj 120071 . l2008al ). these Zeeman signatures fea- 
ture a markedly non-antisymmetric shape (with a blue lobe 
stronger and narrower than the red one) indicating that the 
line is formed within a region featuring a strong velocity gra- 
dient, nicely supporting the idea that it traces non-static, 
chromospheric post-shock plasma decelerating towards the 



surface of the star. Note that longitudinal fields at similar 
phases but different rotation cycles grossly agree with each 
other, even when He I emission fluxes (and thus presum- 
ably accretion rates) are significantly different (e.g., at cycles 
1.609 and 2.542), indirectly confirming that He I emission is 
probing localised regions with well defined field strengths 
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and orientations at the surface of the star. These longitu- 
dinal field estimates are signi ficantly larger than those esti- 
mated in the previous study ijDonati et al.ll20 0~7r). by a typ- 
ical factor of 2 — 3; it reveals in particular that the large- 
scale field of V2f29 Oph strongly varied on a timescale of 
only 4 yr. Maximum longitudinal fields from LSD and He I 
profiles are found to occur at roughly the same phase in the 
present run, while they were shifted by about 0.25 cycle in 
our previous data set; this is additional confirmation that 
the large-scale field topology significantly changed between 
the two epochs. 

Core emission in Ca II 1RT lines is another very useful 
probe of magnetospheric accretion, even though it often fea- 
tures a significant contribution from non-accreting chromo- 
spheres and is thus a more ambiguous proxy than He I emis- 
sion. However, the redder spectral location, the higher mag- 
netic sensitivity and the multiple nature of the correspond- 
ing spectral lines more than compensate for this drawback; 
with their classical, nearly-antisymmetric profiles (with re- 
spect to the line center, see Fig. [1] bottom right panel), the 
Zeeman signatures of Ca II IRT lines are also simpler and 
easier to model than those of the He I line. To extract the 
core Ca n emission profiles (used later in the imaging analy- 
sis, see Sec. O from our spectra, we start by constructing a 
LSD-like weighted average of the 3 IRT lines; we then sub- 
tract the underlying (much wider) photospheric absorption 
profile, with a Lorentzian fit to the far wings over a veloc- 
ity interval of ±200 kms" 1 about the emission core. IRT 
emission is centred at a RV of about —6 kms -1 , i.e., red- 
shifted by only 1 kms -1 with respect to LSD photospheric 
lines (as opposed to 7 kms -1 for the He I line), confirm- 
ing that it traces slowly-moving regions of the post-shock 
accretion funnels that are closer to the surface of the star. 
In 2009 July, equivalent widths of the Ca II emission are in 
average equal to about 18 kms -1 (i.e., 0.050 nm) and vary 
from 12 to 24 kms" 1 (i.e., 0.035 - 0.070 nm) as a result 
of both rotational modulation and intrinsic variability (see 
Fig- right column third panel); maximum IRT emission is 
found to occur at phase 0.70, i.e., 0.15 rotation cycle later 
than He I emission and magnetic maxima. RVs of IRT emis- 
sion are found to vary by about ±2 kms -1 about the mean, 
consistent with a localised emission region centred at phase 
0.70 and at high latitudes. 

Clear Zeeman signatures are detected at all epochs in 
conjunction with Ca n IRT emission, featuring peak-to-peak 
amplitudes of up to about 10% (see Fig. [T] bottom right 
panel). Corresponding longitudinal fields vary from 0.4 to 
1.2 kG (see Fig. [2] left column third panel); this is about 
twice larger than in our previous data set (|Donati et al.l 
2007), confirming again that the field has considerably 
strengthened within the corresponding 4 yr. Maximum lon- 
gitudinal field is observed at phase 0.55, i.e., in conjunction 
with the magnetic maximum traced by He I lines. The longi- 
tudinal field curve is found to repeat fairly well over the two 
rotation cycles, clearly demonstrating that rotational mod- 
ulation is dominant. The ratio of He I to Ca n longitudinal 
fields is ~ 4 across the whole run; this essentially reflects 
that Ca n longitudinal fields suffer from dilution and flux 
cancellation (with respect to those from He i) as a result of 
Ca II emission forming over a much wider area (including 
in particular the quiet chromosphere) while He I emission 
traces only the accretion regions. 



As usual for cTTSs, Balmer lines (and in particular Ha 
and H/3) are in emission and show clear temporal variability 
and rotational modulation (see Fig.01ower panel). As in our 
previous data set, Ha and H,9 mostly show one main emis- 
sion peak centred at about -15 kms -1 , blue-shifted by about 
8 kms -1 with respect to the LSD photospheric profiles; the 
average equivalent widths of both profiles are equal to 800 
and 150 kms -1 respectively (i.e., 1.7 and 0.24 nm), i.e., 
very similar to those measured in 2005 June. Looking at the 
variance profiles of both lines (see Fig. [4j, one can see that 
temporal variations mostly affect the central emission peak 
and the red wing, both parts of the line being anti-correlated 
with each other; this is similar to what was reported from 
our previous data set, except that variations in the central 
peak are now smaller than in 2005 June while those in the 
red wing are both larger and shifted to smaller velocities 
(with the variance profile now peaking at 100 kms -1 in- 
stead of 150 kms" 1 ). In our new data set, Balmer emission 
and profile variations are also clearly visible in the blue wing 
(around —170 kms -1 ), both in Ha and H/3 (see Fig. [4]); this 
emission component was not detected in the previous data 
set. 

Recurrent absorption episodes in the red wing of Balmer 
lines are usually interpreted as evidence for accretion veils 
anchored in the inner regions of the surrounding disc and 
perio dically intersecting the line of sight as the star rotates 
fe.g.. lBouvier et al.ll2007l l. Those observed on V2129 Oph in 
2009 July occur between phases 0.6 and 0.9 (see Fig.[3]lower 
panel), i.e., 0.05 — 0.35 rotation cycle (or 0.20 cycle in aver- 
age) after He I and magnetic maxima; it suggests that the 
corresponding accretion veils are trailing the main magnetic 
spot detected at the surface of the protostar (centred at 
phase 0.55). In 2005 June, these absorption episodes were 
shorter, being visible for less than 15% of a rotation cy- 
cle and occurring sli ghtly (about 0.05 ro tation cycle) before 
magnetic maximum jPonati et alj|2007l ). 

From the average equivalents widths of the He I, Ca II 
IRT, H/3 and Ha lines, we derive logarithmic line fluxes (with 
respect to the luminosity of the Sun L©) equal to —5.1, 
—4.8, —4.1 and —3.2 respectively, implying logarithmic ac- 
cretion luminosities (with respect to Lq) of —2.0, —2.0, —2.3 
and —1.7 respe ctively (using empirical correlations from 
iFang et al.ll2009h . As estimates from the two main accretion 
proxies agree with each other and with the overall mean, 
we can safely conclude that the average logarithmic accre- 
tion luminosity of V2129 Oph is —2.0 ± 0.3 and thus that 
the average logarithmic mass accretion rate (in Moyr -1 ) 
is equal to —9.2 ± 0.3. Using the same relations, we find 
that the logarithmic accretion rate at magnetic maximum is 
typically —9.0 (with peaks of up to —8.9) and about —9.4 
around magnetic minimum (i.e., phase 0.05), the observed 
fluctuation most likely tracing the varying viewing config- 
uration (rather than an intrinsic change in the accretion 
rate). Mass accretion rates can independently be estimated 
through the full width of Ha at 10% height fe.g.. lNatta et al.l 



2 To derive line fluxes from normalized equivalent widths, we ap- 
proximate the continuum level by a Planck function at the tem- 
perature of the stellar photosphere. Results are fo und to be com- 
patib le with those in the published literature (e.g. jMohantv et al.l 
120051) within better than 0.1 dex. 
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Figure 3. Variation of the He I D3 (Stokes /: top left, Stokes V: top right), Ha (bottom left) and H/3 (bottom right) lines of V2129 Oph 
in 2009 July. To emphasize variability, the average profile over the run is shown in red. Rotation cycles (as listed in Table 1) and 3<r 
error bars (for Stokes V data only) are also mentioned next to each profile. 
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Ha (left panel) and H/3 (right panel) of V2129 Oph at both observing epochs. Standard deviation profiles are expanded by a factor of 3 
and shifted upwards by 1 for graphics purposes. 



12004 ICieza et aUl201Ch ; in the case of V2129 Oph in 2009 
July, Ha exhibits a full width of 350 kms -1 implying an 
logarithmic mass accretion rate estimate of —9.4 ± 0.6 in 
reasonably good agreement with our main estimate. 

Line equivalent widths and fluxes measured during our 
2005 run (and in particular those of He I and Ca n IRT ac- 
cretion proxies) are basically identical to those listed here, 
suggesting that the mass accretion rate at the surface of 
V2129 Oph is roughly stable (within a factor of ~2) on long 
timescales as well; in particular, this is different than the 
behaviour reported for AA Tau, where the mass accretion 
rate can vary b y an order of magni tude between successive 
rotation cycles |Donati et al.ll2010bl ). Further data collected 
over a longer time span are however needed to confirm that 
we did not catch V2129 Oph in an unusually quiet state of 
accretion at both epochs. 



emission map compatible with the series of rotationally mod- 
ulated Stokes I and V LSD and Ca n IRT profiles. The re- 
construction process is iterative and proceeds by comparing 
at each step the synthetic Stokes I and V profiles corre- 
sponding to the current images with those of the observed 
data set. The magnetic field is described through its poloidal 
and toroidal com ponents expressed as spherical-harmonics 
(SH) expansions l|Donati et al.l [20061 ). The spatial distribu- 
tion of photospheric brightness (with respect to the quiet 
photosphere) and that of accretion-powered Ca n emission 
(in excess of and with respect to that produced by the quiet 
chromosphere) are modelled as series of independent pixels 
(typically a few thousand) on a grid covering the visible sur- 
face of the star (with spots in the brightness image assumed 
to be darker/cooler than the quiet photosphere and spots 
in the accretion-powered Ca n emission map supposed to be 
brighter than the quiet chromosphere). 



5 MODELLING THE SURFACE OF V2129 OPH 

5.1 Overview of the modelling method 

Using the set of LSD and Ca n IRT profiles of V2129 Oph 
described in the previous section, we can recover the large- 
scale field topology at the surface of the protostar, as well as 
maps of how photospheric brightness and accretion-powered 
Ca n emission distributes with longitude and latitude. In 
this aim, we apply our new modelling technique, detailed 
extensively and prov ed successful in a previous MaPP study 
(|Donati et al.ll2010bl ). In the present paper, we only briefly 
describe the technical aspects of this method and concen- 
trate on its specific application to our new V2129 Oph data 
set; we also carry out a new analysis of the older data set 
(collected in 2005 June and a nalysed with an ex perimen- 
tal version of our imaging tool, iDonati et aT]|2007l ) to allow 
a more direct and meaningful comparison of reconstructed 
maps. 

Following the principles of maximum entropy, our code 
automatically derives the simplest magnetic topology, pho- 
tospheric brightness image and accretion-powered Ca II 



Synthetic profiles are computed by summing up the el- 
ementary spectral contributions from each image pixel over 
the visible hemisphere, taking into account all relevant local 
parameters of the corresponding grid cell (e.g., brightness, 
accretion-powered excess emission, magnetic field strength 
and orientation, radial velocity, limb angle, projected area). 
Since the problem is partly ill-posed, we stabilise the inver- 
sion process by using an entropy criterion (applied to the 
SH coefficients and to the brightness/excess emission image 
pixels) aimed at selecting the field topologies and images 
with minimum information among all those compatible with 
the data. The relative weights attributed to the various SH 
modes can be imposed, e.g., for purposedly producing an- 
tisymmetric or symmetric field topologies with respect to 
the centre of the star (by favouring odd or even SH modes, 
IDonati et ai1l2007l . l2008al ) . More details concerning the spe- 
cific des cription of local profi les used in the model can be 
found in IDonati et all |2010bl ). 
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Figure 5. Measured (red open circles) and fitted (blue pluses) equivalent widths of the Ca II IRT LSD profiles of V2f29 Oph in 2009 
July (left panel) and 2005 June (right panel). The model wave (green line) providing the best (sine+cosine) fit to the data (with a period 
of 6.53 d) presumably traces rotational modulation, while the deviation from the fit illustrates the strength of intrinsic variability. 



5.2 Application to V2129 Oph 

Our imaging model assumes that the observed profile vari- 
ations are mainly due to rotational modulation, and po- 
tentially to surface differential rotation as well when the 
star is observed for at least several rotation cycles; all other 
sources of profile variability (and in particular intrinsic vari- 
ability) cannot be properly reproduced and thus contribute 
as noise into the modelling process, degrading the imaging 
performance and potentially even drowning all relevant in- 
formation. Filtering out significant intrinsic variability from 
the observed profiles is thus worthwhile to optimise the be- 
haviour and convergence of the imaging code. 

We implement this by applying specific corrections on 
our data set. We first suppress veiling by scaling all LSD 
Stokes / and V photospheric profiles, to ensure that un- 
polarized lines have the same equivalent widths. We also 
remove the non rotationally-modulated part in the observed 
fluctuations of Ca II IRT emission, by fitting them with a 
sine+cosine wave (see Fig. [5} and by scaling the correspond- 
ing Stokes / and V profiles accordingly, thus ensuring that 
equivalent widths of unpolarized lines match the optimal fit. 
Although only approximate (especially the removal of the 
intrinsic variability) , this procedure is at least very straight- 
forwa rd and has proved su ccessful when applied to real data 
(e.g., iDonati et ail l2010bl ) and efficient at retaining rota- 
tional modulation mostly. 

The reconstructed magnetic, brightness and accretion 
maps of V2129 Oph are shown in Fig. [5] with correspond- 
ing fits to the data shown in Figs. [7] and [8] we also include 
new maps (and fits) derived from the previous data set, but 
reconstructed with our new modelling method to ease com- 
parison between the two sets of images. The SH expansion 
describing the field was limited to terms with I < 8, very 
small power being reconstructed in higher order modes. Both 
images were reconstructed assuming that the field is anti- 
symmetric with respect to the centre of the star (see below) . 

The overall fits to the data are reasonably good, repro- 
ducing the data down to reduced chi-squares Xi of about 2 
(starting from initial Xi of 25 and 40 for the 2009 July and 
2005 June data sets respectively) . Optimal fits are obtained 



for v sin i = 14.5 ± 0.5 kms -1 and « ra d = — 7.0 ± 0.1 kms -1 ; 
Ca II IRT lines are found to be slightly redshifted (by about 
1 kms -1 ) with respect to photospheric lines, as usual for 
cTTSs. We also find that optimal fits to the observed Stokes 
V profiles can be achieved when setting the local filling fac- 
tor ip, describing the relative proportion of magnetic areas 
at any given point of the stellar surface, to V — 0.4 ±0.1. 
The emission profile scaling factor e, describing the emis- 
sion enhancement of accretion region s over the quiet chro - 
mosphere, is again set to e = 10 (as in lDonati et al.ll2010tj ). 
yielding fractional areas of accretion spots of about 2 — 3% 
at both epochs. Most other model parameters, and in partic- 
ular those describing the local line profiles, are either iden- 
tical of very similar to those chosen in our pre vious studies; 
for mor e information, the reader is referred to lDonati et al.l 
(|2010bl ). 

5.3 Results 

The reconstructed large-scale magnetic topology of 
V2129 Oph in 2009 July is mostly poloidal, with only 5% of 
the magnetic energy in the toroidal component; this poloidal 
component is mainly antisymmetric, with at least ~ 60% 
of the energy concentrating in modes with m < 1/2. We 
find that the field is more likely to be antisymmetric about 
the centre of the star, a symmetric topology requiring much 
more magnetic energy to produce the same fit to the data. 
The dominant terms in the reconstructed SH expansion of 
the large-scale poloidal field correspond to I = 3 (with 
m — 0, 1, i.e., a moderately tilted octupole) and £ — 1 (with 
m = 0, i.e., an aligned dipole), whose respective intensity at 
the pole are about 2.1±0.3 kG and 0.9±0.1 kG respectively; 
from different reconstructions with slightly different input 
parameters and reconstruction schemes, we estimate that 
the octupole is typically 2 — 3 stronger than the dipole. The 
main magnetic feature that we reconstruct is a high-latitude 
radial field region where the magnetic intensity peaks at 
about 4 kG, slightly smaller though comparable to the max- 
imum fields traced by the He I line (about 5 kG, see Fig. 0; 
this main magnetic region is centred at phase 0.5, compati- 
ble with the estimate guessed in Sec.|4]from the longitudinal 
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Figure 6. Maps of the radial, azimuthal and meridional components of the magnetic field B (first and third rows, left to right panels 
respectively), photospheric brightness and excess Ca II IRT emission (second and fourth rows, first and second panels respectively) at 
the surface of V2129 Oph in 2009 July (two upper rows) and 2005 June (two lower rows). Magnetic fluxes are labelled in G; local 
photospheric brightness (normalized to that of the quiet photosphere) varies from 1 (no spot) to (no light); local excess Call emission 
varies from (no excess emission) to 1 (excess emission covering 100% of the local grid cell, assuming an intrinsic excess emission of 
lOx the quiet chromospheric emission). In all panels, the star is shown in flattened polar projection down to latitudes of —30°, with the 
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Figure 7. Maximum-entropy fit (thin red line) to the observed (thick black line) Stokes / and Stokes V LSD photospheric profiles (first 
two panels) and Ca II IRT profiles (last two panels) of V2129 Oph in 2009 July. The light-blue curve in the bottom left panel shows the 
(constant) contribution of the quiet chromosphere to the Stokes / Ca II profiles. Rotation cycles and 3cr error bars (for Stokes V profiles) 
are also shown next to each profile. 
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Figure 8. Same as Fig.[7]for the 2005 June data set. 



field curves of both Ca II IRT and He I lines. The octupole 
component of the poloidal field can be directly traced in the 
reconstructed magnetic map, as a low-latitude ring of nega- 
tive radial field encircling the main high-latitude radial field 
region. 

The reconstructed magnetic maps of V2129 Oph in 2005 
June are grossly similar to those of 2009 July except for the 
large change in the overall field intensi ty; they are al s o com - 
patible with the preliminary version of fponati et~all \200H ). 
Again, the large-scale field is found to be mostly poloidal 
and axisymmetric, with the poloidal component being dom- 
inated by a moderately tilted octupole (with a polar strength 
of about 1.5 ± 0.3 kG); the dipole component is typically 
5 times smaller than the octupole component, with polar 
strengths of about 0.3 kG. The high-latitude radial field re- 
gion features peak intensities of about 2 kG, slightly larger 
though compatible with fields traced by the He I line. The 
phase at which this high field region is located is different 
in 2009 July and 2005 June; however, this difference likely 
reflects no more than the uncertainty on the rotation pe- 
riod. The reconstructed maps thus confirm the preliminary 
conclusions of Sec. [4] about the overall strengthening of the 
field (by typically a factor of 2) between 2005 and 2009; they 
further suggest that this increase mostly concerns the dipole 
component of the field (about thrice stronger in 2009) and to 
a lesser extent its octupole component (about 50% stronger 
in 2009). 

The photospheric brightness map we reconstruct for 
V2129 Oph in 2009 July essentially features one large dark 
spot at high latitude covering about 6.5% of the total stellar 
surface (i.e., about 25% of the visible hemisphere at maxi- 



mum visibility) and within most of which the brightness con- 
trast at visible wavelengths (with respect to the surround- 
ing photosphere) is larger than 20; this dark spot mostly 
overlaps the strongly magnetic region discussed above, mak- 
ing the latter hardly visible at optical wavelengths (hence 
its non-detection in LSD profiles and in the correspond- 
ing longitudinal field curve, see Fig. [2}. This photospheric 
brightness map is also grossly consistent (both in phase 
and amplitude) with the observed photometric light-curve 
of V2129 Oph, even though the latter was not used in the 
modelling. The photospheric brightness map corresponding 
to 2005 June (again compati ble with, though le ss detailed 
than, the previous version of iDonati et al.ll2007t ) features a 
smaller and less contrasted high-latitude dark spot as well 
as low-latitude appendages, covering altogether about 5% of 
the total stellar surface. As in 2009 July, the high-latitude 
spot is mostly coincident with the strongly magnetic region 
of the radial field map. We estimate that the increase in size 
and contrast of the high-latitude dark spot on V2129 Oph 
between 2005 and 2009 is real, and readily traceable in the 
comparatively larger distortions (by about a factor of 2 in 
average, as judged from standard deviation profiles) in the 
unpolarized LSD profiles of the 2009 data set (see first pan- 
els of Figs. [7] and 0). 

The map of accretion-powered excess emission that we 
reconstruct from the 2009 July data set shows one main 
high-latitude region, covering about 2.5% of the total stellar 
surface. This region is located at phase 0.6, i.e., slightly trail- 
ing by about 0.1 cycle the main radial field region and the 
dark photospheric spot; this phase shift is readily traceable 
in the rotational modulation of both the equivalent widths 
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and radial velocities of the Ca II IRT emission (see Sec. |4}, 
suggesting that it is real even though corresponding to a 
spatial size comparable to our resolution element at the sur- 
face of the star. A similar accretion region is reconstructed 
from the 2005 June data set, this time overlapping the main 
field region and dark photospheric spot. 

Our modelling can also estimate the amount by which 
the photosphere i s shear ed by differential rotation (e.g., 
iDonati et al"1l2003l . l2010al ). This is done by assuming that 
the rotation rate at the surface of the star is varying with 
latitude 8 (as Q cq — dQ, sin 2 8 where f2 oq is the rotation rate 
at the equator dQ the difference in rotation rate between the 
equator and the pole), and by determining the pair of dif- 
ferential parameters Q eq and dtt that produces the best fit 
to the data at a given image information content. Since dif- 
ferential rotation is estimated from subtle departures with 
respect to pure rotational modulation, we only use LSD pho- 
tospheric profiles for this task, in order to minimise all other 
potential sources of intrinsic variability. When using Stokes 
V profiles, we find a clear minimum in the x? map, located at 
fi cq = 0.974 ±0.005 radd -1 and dQ. = 0.036 ±0.015 radd" 1 
(see Fig. [9J. The corresponding rotation periods for the 
equator and pole are 6.45 ± 0.04 d and 6.70 ±0.15 d re- 
spectively, compatible (within ~ 2cr) with the range of re- 
ported photometric periods (ranging from 6.35 to 6.60 d 
iGrankin et al.| [2008V This is also compatible (within ~ la) 
with the preliminary estimate derived from the 2005 Stokes 
/ data set. No clear x? minimum is observed when using the 
LSD Stokes J profiles from 2009 July. 

From our new estimate alone, we cannot firmly exclude 
that V2129 Oph is rotating rigidly; if this were the case, the 
most likely rotation period (i.e., that minimising x?) would 
be 6.52 ±0.01 d, close to that assumed here for phasing data 
(equal to 6.53 d, see Eq. fJJ). Rigid rotation is however un- 
likely given the observed variability of the photometric pe- 
riod. At 99% confidence level, we can ascertain that the pho- 
tospheric shear of V2129 Oph is Sun-like in sign (i.e. with a 
fast equator and slow pole) and smaller than ~ 0.07 radd -1 
(about that of the Sun). The corresponding time for the 
equator to lap the pole by one complete rotation cycle is 
found to be 175^52 5 d, and is larger than 90 d at 99% con- 
fidence level. Images presented in Fig. [15] are reconstructed 
assuming our new estimate of the photospheric shear; im- 
ages derived assuming rigid rotation (with a period of 6.53 d) 
are however virtually identical to those presented here (as 
expected from the small differential rotation). 



6 SUMMARY AND DISCUSSION 

This paper presents new results of the MaPP project, fo- 
cussing on the impact of magnetic fields on star formation, 
during the cTTS phase in particular. Here we concentrate on 
V2129 Oph, the brightest cTTS of the p Oph star formation 
cloud, for which we present new spectropolarimetric data 
collected in 2009 July with ESPaDOnS@CFHT, covering 2 
complete rotation cycles of the protostar and augmented 
from almost simultaneous photometric data collected with 
SMARTS; additional multiwavelength data (e.g., X-ray data 
from Chandra/HETGS, visible spectra from HARPS@ESO 
and SMARTS) were also collected as part of this campaign 
and will be presented in forthcoming companion papers. 
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Figure 9. Variation of x? as a function of f2 eq and dQ, derived 
from the modelling of Stokes V LSD profiles at constant informa- 
tion content. The outer color plot corresponds to a x 2 increase 
of 1% corresponding to a 99% confidence interval (for each pa- 
rameter considered separately) given the number of data points 
adjusted in the process (644). 



From these time-resolved data and using the last version 
of our tomographic techniques, we reconstruct maps of the 
large-scale magnetic field, of the photospheric brightness 
and of the accretion-powered Ca II IRT excess emission at 
the surface of V2129 Oph; we also include in the paper a 
complete re-analysis of our old spectropolarimetric data col- 
lected in 2005 July, to ease comparison with our newest re- 
sults. 

The large-scale field of V2129 Oph is found to be mostly 
poloidal and axisymmetric, and essentially consists of a 
~2.1 kG octupole and a ~0.9 kG dipole, both tilted by about 
20° with respect to the rotation axis; the main magnetic fea- 
ture reconstructed on the map is a 4 kG high-latitude radial 
field region. We also find that the large-scale field strongly 
increased by typically a factor of 2 since 2005 June, with 
the dipole component enhanced by a factor of ~3 and the 
octupole component by about 50%; in this process, the oc- 
tupole to dipole intensity ratio has decreased by a factor 
of ~ 2, from about 5 (in 2005) down to ~ 2.5 (in 2009). 
This change is readily visible on the main high-latitude ra- 
dial field feature, whose average magnetic strength (as de- 
termined from, e.g., Ca II IRT or He I lines) has increased by 
at least a factor of 2 between both epochs. We finally report 
that the surface differential rotation by which the detected 
large-scale field is sheared is Sun-like in sign (with a fast 
equator and slow poles) and weaker than solar in strength 
(by about a factor of 2). 

We consider that the change in the large-scale magnetic 
topology of V2129 Oph is unambiguous evidence that its 
field is generated by dynamo processes, rather than being 
a fossil remnant of earlier formation stages; although this 
conclusion was anticip ated both on observational and the- 
oretical grounds (e.g., IDonati fc Landst rect 2009]), its firm 
demonstration is nevertheless a significant step forward in 
our understanding of the origin of magnetic fields in low- 
mass protostars. Our result suggests that dynamo fields of 
cTTSs are variable on a timescale of a few years and are thus 
intrinsically non-stationary; a regular magnetic monitoring 
of a small sample of low-mass protostars should be able to re- 
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veal whether these dynamos are cyclic, and if so, provide an 
estimate of the typical period of their magnetic cycles. The 
recent report that two cTTSs with virtually identical stel- 
lar parameters and evolutionary states (namely BP Tau and 
AA Tau) show discrepant magnetic topologies (with unequal 
dipole to octupole intensity ratios, respectively equal to ~1 
and >5 respectively) and different dipole strengths in par- 
ticular (by about a factor of 2), although intriguing at first 
sight, seems consistent with what we observe on V2129 Oph 
(where we see a comparable change in dipole field strength 
between 2005 and 2009). 

Our data also indicate that areas of enhanced accretion- 
powered Ca n and He I emission are close to the main mag- 
netic region of V2129 Oph, as in all cTTSs monitored with 
MaPP up to now; it confirms in particular that this strong 
field region (and its immediate vicinity) is where accretion 
funnels, linking the protostar to the inner edge of its ac- 
cretion disc, are anchored. We further demonstrate that in 
V2129 Oph, these strong magnetic regions are coincident, 
at photospheric level, with large dark spots featuring a high 
brightness contrast with the neighbouring photosphere and 
generating most of the observed photometric variations. We 
also find that this dark spot is increasing in size and con- 
trast between 2005 and 2009, i.e., when the field in main 
magnetic region strengthens, in line with theoretical expec- 
tations that stronger fields are more efficient at inhibiting 
local convection. We note that the observed photometric 
modulation of V2129 Oph cannot be attributed to partial 
occultations of the protostar by the accretion warp at the in- 
ner d isc edge like those seen on AA Tau (e.g.,H ouvier et al.l 
2007); these occultations are indeed generating photomet- 
ric eclip ses that are deep (u sually larger than f mag), color 
specific feouvier et al.|[l999T) and lagging phases of magnetic 
maxima l|Donati et aTE oiObl. at variance with what we re- 
port here for V2f 29 Oph. It confirms at the same time that 
the rotation axis of V2129 Oph is likely not inclined by more 
that 60° with respect to the line of sight (see Sec.[2|. 

Using the He I D3, Can IRT, H/3 and H a emission 
fluxes along with the empirical correlations of iFang et al.l 
(2009), we derive that the average logarithmic accretion 
rate of V2I29 Oph (in M yr _1 ) is -9.2 ± 0.3, peaking at 
—9.0 in average at magnetic maximum (i.e., phase 0.55). 
It suggests that the mass accretion rate of V2129 Oph is 
stablqj to within a factor of ~2 on both short and long 
timescales (i.e., weeks to years). Given the intensity of the 
dipole component reconstructed for V2129 Oph (equal to 
B„ ~ 450 G and ~ 150 G at the equator in 2009 and 2005 
respectively, see Sec. [5]) and using a logarithmic mass accre- 
tion rate of —9.0, we obtain that the radius r mag at which the 
inner disc of V2129 Oph is truncated is r mag ~ 7.2 and 
~ 3.9 -R* for 2009 and 2005 respectiveljQ, or equivalently 
that r mag /r cor ~ 0.93 and ~ 0.50 at the same two e pochs 
(using the theoretical estimates of lBessolaz et al.ll2008h . 



3 More data collected over a longer time span are however needed 
to confirm that we did not catch V2129 Oph in an unusually quiet 
state of accretion at both epochs. 

4 Note that our previous estimate of r mag in 2005 June 
jDonati et al.|[2007l ). derived from a different expression of r mag 
and a different mass accretion rate, should be considered as ob- 
solete and superseded by the present one. 



This shows that, in 2009 July, V2129 Oph has a mag- 
netosphere extending out to r mag ~ r COI , where the dipole 
component of the field vastly dominates that of the octupole 
(by a factor of about 25 in field strength) even though the 
latter dominates the former at the stellar surface (by about 
a factor of 2.5 and 1.9 at the poles and equator respec- 
tively); in this situation, accretion is expected to occur nat- 
urally towards the high-latitude magnetic poles only (and 
not towards low- latitudes, as in a pure octupole case, e.g., 
lLong et al l 120091 1. Our observations confirm that this is in- 
deed the case, with mass accretion at the surface of the star 
concentrating at a magnetic colatitude smaller than 10° (see 
Fig- O ■ I n 2005 June, we find that the magnetosphere only 
extends to r mag ~ 3.9 R+ while at the same time the dipole 
field is half as strong as in 2009 with respect to the octupole 
component (the octupole dominating the dipole by a factor 
of 5 at the poles and 3.8 at the equator, at the surface of the 
star); we nevertheless find that the dipole still dominates 
the octupole at r mag , but only by a factor of about 4 in field 
strength. We speculate that this is enough for surface accre- 
tion to concentrate mostly at magnetic poles only (rather 
than at low latitudes) as suggested again by observations. 

Our results are in relatively good agreement with those 
of recent very detailed 3D MHD simulations carried ou t in 
the specific case of V2129 Oph (|Romanova et alj|20ich . In 
particular, estimates of the relative extent of r mag with re- 
spect to r cor and its dependence with dipole strength and 
accretion rate are roughly compatible. The grossly circu- 
lar shape and polar location of accretion spots found in 
our maps is also similar to theoretical predictions when 
r mag /r C or > 0.7 and differ in particular from the pre- 
dicted crescent-shape d accretion regions in purely dipolar 
magnetospheres fe.g.. iRomanova et al.ll2003l ); the accretion 
arcs/rings that simulations predict at lower latitudes when 



? /r c 



0.5 (i.e., as in 2005) are however not observed, 



for a yet unclear reason. 

Balmer profiles further document the different magneto- 
spheric accretion configurations of V2129 Oph in 2009 July 
and 2005 June. The standard deviation profiles clearly show 
that the red-shifted absorption in 2009 July occurs at sig- 
nificantly smaller velocities than in 2005 June (see Fig. Q. 
The origin of this change is not clear yet; multiple observa- 
tions are needed to investigate this issue in more details. In 
addition, we observe that this red-shifted Balmer absorption 
is trailing magnetic maximum in 2009 July (by about 0.15 
cycle in average, see Sec. 2]), while the opposite is observed 
in 2005 June. We suggest it indicates that part of the accre- 
tion funnel is anchored beyond r cor in 2009 July (as expected 
from r mag ~ r cor ) with the corresponding flux tubes lagging 
the main funnel as they are dragged behind by the slow Kep- 
lerian flow; since r mag <C r cor in 2005 June, we may logically 
expect an opposite situation (with the base of accretion fun- 
nels being dragged ahead of the main magnetic region by the 
fast Keplerian flow), in rough qualitative agreement with ob- 
servations. Finally, blue-shifted Balmer emission is detected 
at velocities of ~ —170 kms -1 in 2009 July while none was 
seen in 2005 June, suggesting that V2129 Oph is now trig- 
gering a significant outflow; we propose that this outflow 
comes from the inner disc, tracing the (small) fraction of 
the accreting disc material expelled outwards, as expected 
when r mag ~ r cor . 

We also report that the high-latitude accretion-powered 
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Ca II IRT excess emission spot that we detect at the surface 
of V2129 Oph (at phase 0.6 - 0.7) is slightly trailing the 
main magnetic region as well as the accretion-powered He I 
emission spot (both centred at phase 0.5) by about 0.1 cy- 
cle in 2009 July (see Sec. Hf, while all three regions were 
more or less overlapping in 2005 June. The origin of this 
small (but apparently real) spatial shift is not clear yet. We 
suspect that this phenomenon relates to the trailing accre- 
tion veils reported above, even though red-shifted Balmer 
absorption obviously traces a physically different magneto- 
spheric plasma than Ca II IRT excess emission; we propose 
for instance that it reflects the complex temperature struc- 
ture of the postshock region that results from the sheared 
accretion funnels. More data are obviously needed to inves- 
tigate this issue in a quantitative way. 

As a summarising conclusion, we propose a schematic 
description of how magnetic fields are likely to affect the evo- 
lution of cTTSs, and in particular their angular momentum. 
Our new spectropolarimetric data set on V2129 Oph shows 
that magnetic fields of cTTSs exhibit a large-scale struc- 
ture, and in particular a dipolar component, that varies on a 
timescale of a few years, as a likely result of a non-stationary 
dynamo. In this context, magnetospheric accretion, essen- 
tially controlled by the lowest orders of the large-scale field 
of the protostar, is expected to reflect more or less directly 
changes in the magnetic topology. When the dipole compo- 
nent of the field is strong enough to ensure that r mag > r cor , 
the protostar enters a propeller-like regime and is actively 
spun-down; when the dipole field gets weaker, r mag shrinks 
back within r COI and the star spins up again. Since the typi- 
cal fluctuation timescale of dynamo fields (of order of years 
to decades) is much shorter than that of star-disc angu- 
lar momentum transfer (by at least 4 orders of magnitude), 
the protostar ends up rotating at an average rate; this rate 
more or less reflects the mean dipole strength (with larger 
fields implying slower rotation) and mean mass accretion 
rate (with stronger accretion implying faster rotation) over 
the averaging timescale. As a result, cTTSs with masses in 
the range 0.5 — 1.3 Mq rotate slowly, their dynamo fields 
being stro ngly dipolar and ax isymmetric like those of mid- 
M dwarfs (|Morin et al.ll2008l) ; cTTSs with either lower or 
larger masses rotate more rapidly, being appa rently less ef- 
ficient at producing strong dipole fields (e.g., iDonati et al.l 
1201031 ). A A Tau and BP Tau obviously belong to the first 
catego ry, both host i ng a dipole c omponent stronger than 
1 kG (Donati et all l2008al . l2010bh . With a radiative core 
growing in size rapidly, V2129 Oph is apparently no longer 
successful at building up a strong dipole field (or no more 
than for a small amount of time) and can likely no longer 
counteract its spin-up. 

Upon completion, MaPP should be able to assess (on 
a small statistical sample) how well our proposed picture 
matches cTTSs; in particular, repeated and regular visits 
on a limited, well-selected subgroup, like that presented here 
for V2129 Oph, will be essential to address the important 
issue of non-stationary dynamos in low-mass protostars, to 
evaluate the constant and fluctuating parts of their dipolar 
fields, and to quantify the associated impact on their angular 
momentum evolution. 
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